Chapter 5 

Broken Symmetry and Collective 
Phenomena 

Previously, we have seen how the field integral method can be deployed in many-particle 
theories. In the following chapter, we will learn how elements of perturbation theory can be 
formulated efficiently by staying firmly within the framework of the field integral. In doing 
so, we will see how the field integral provides a method for identifying and exploring non- 
trivial reference ground states — 'mean-fields '. A fusion of perturbative and mean-field 
methods will provide us with analytical machinery powerful enough to address a spectrum 
of rich applications. 

Historically, the effects of interaction on many-body systems are typically dealt with 
within the framework of "diagrammatic" perturbation theory, a series expansion in the 
interaction strength. However, in the following, we will search for a different methodology. 
Our motivations are two-fold: Firstly, the structures that typically appear from pertur- 
bative expansions can be assimilated more straightforwardly. But, more importantly, the 
development of phase instabilities typically reflect the appearance of non-perturbative 
structures. Thus, what we would like to develop is a theoretical framework that is ca- 
pable of (a) detecting the 'right' reference states or 'mean-fields' of a system, (b) that 
enables us to efficiently apply perturbative methods around these states and, finally, (c) 
to do this in a language that draws upon the 'physical' rather than the plain microscopic 
degrees of freedom as the fundamental units. 

To this end, in the following sections we will develop a functional integral based ap- 
proach that meets these criteria. In contrast to the previous chapters, the discussion here 
will be decidedly biased towards concrete application to physically motivated problems. 
After the formulation of the general strategy of field integral based mean-field methods, 
the next section will address the problem of the weakly interacting electron gas. The 
exemplification of the new concepts on "familiar territory" will enable us to better un- 
derstand the intimate connection between the mean-field approach and straightforward 
perturbation theory. In subsequent sections we will then turn to the discussion of problems 
which lie firmly outside the range of direct perturbative summation, viz. the phenomena 
of superfluidity and superconductivity. 
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5.1 Mean-Field Theory 



Roughly speaking, the functional approach to problems with a large parameter proceeds 
according to the following programme: 

1. In the first place, one must identify the relevant structural units of the theory. (This 
part of the programme can be efficiently carried out by the straightforward methods 
discussed earlier.) 

2. Secondly, it is necessary to introduce a new field — let us call it <fi for concreteness 

- that encapsulates the relevant degrees of freedom of the low energy theory. 

3. With this in hand, one can then trade integration over the 'microscopic fields' for 
an integration over 0, a step often effected by an operation known as the Hubbard- 
Stratonovich transformation. 

4. The low-energy content of the theory can often be explored by subjecting the re- 
sulting action S[(f>] to a stationary phase analysis. (The justification for applying 
stationary phase methods is provided by the existence of a large parameter N >> 1.) 
Often, at this stage, instabilities in the theory show up — an indication of a physi- 
cally interesting problem! 

5. Finally, the nature of the elementary (collective) excitations above the ground state 
can by explored by expanding the functional integral around the solution of the 
stationary phase equations — the 'mean-field'. From this low-energy effective action, 
one can compute physical observables. 

In the next section, we will illustrate how such a programme can be implemented on a 
specific example which can also be studied by perturbative means: 



5.2 Plasma Theory of the Interacting Electron Gas 



In first quantised form, the weakly interacting electron gas is described by the many-body 
Hamiltonian, 



For simplicity, we have chosen to neglect to presence of any underlying lattice potential. 
Formally, the corresponding quantum partition function is obtained as Z = tre~ /3 ( // ~ AtAr - ) 
where the trace runs over a complete basis of many-body states. 

D> INFO. Before plunging into the computation of the quantum partition function from the 
coherent state path integral we should first try to understand in what limit Coulomb interactions 
may be thought of as weak. In fact, the control parameter is the average density of the electron 
gas: n= 1/vq where tq denotes the average interelectron spacing. Coulomb interaction is weak 



N 
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(5.1) 
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when the average potential energy (measured in units of e 2 /ro) is small as compared to the 
typical kinetic energy (measured in units of h 2 /mr^). The ratio of energy scales defines the 
dimensionless density parameter 

e 2 mr 2 . r$ 

— r~ = — = r s, 

where ao = h/e 2 m denotes the Bohr radius. Physically, r s is the radius of the spherical volume 
containing one electron on average; the denser the electron gas, the smaller r s . 

Below, we will be concerned with the limit of high density r s <C 1, in which the effects of 
Coulomb interaction can be treated perturbatively. In the opposite limit, r s S> 1, properties 
become increasingly dominated by the Coulomb interaction. Ultimately, for sufficiently large r s 
(or low density) it is believed that the electron gas undergoes a (first order) phase transition 
to a condensed or 'solid' phase known as a Wigner crystal. (Indeed, this phenomenon is the 
continuum counterpart of the Mott-Hubbard transition described in section 2.2.3.) Although 
Wigner 1 crystals have never been unambiguously observed, several experiments performed on 
low density electron gases are consistent with a Wigner crystal ground state. Quantum Monte- 
Carlo simulation suggests that Wigner crystallisation may occur for densities r s > 37. (Note 
that this scenario relies crucially on being at low temperatures, and the long-range nature of the 
Coulomb interaction. In particular, if the Coulomb interaction is screened V(r) ~ e" r / A , r s ~ 
(?"o/ a o) e_r0// ^ an d the influence of Coulomb interaction at high densities becomes diminished.) 

For r s ~ 0(1), the potential and kinetic energies are comparable. This regime of intermediate 
coupling is notoriously difficult to describe quantitatively. Yet most metals lie in a regime of 
intermediate coupling 2 < r s < 6. Fortunately, there is overwhelming evidence to suggest that a 
weak coupling description holds even well outside the regime over which microscopic theory can 
be justified. The phenomenology of the intermediate coupling regime is the realm of Landau's 
Fermi Liquid Theory. 




Figure 5.1: Sketch showing the adiabatic continuity of the eigenstates in a one-dimensional 
potential well. 

The fundamental principle underlying the Fermi liquid theory is one of "adiabatic conti- 
nuity" [5]: In the absence of an electronic phase transition (such as Wigner crystallisation or 
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the Mott transition), a non-interacting ground state evolves smoothly or adiabatically into the 
interacting ground state as the strength of interaction is increased. 2 An elementary excitation 
of the non-interacting system represents an "approximate excitation" of the interacting system 
(i.e. the 'lifetime' of an elementary excitation is long). Excitations are quasi-particles (and 
quasi-holes) above a sharply defined Fermi surface. The remarkable success (as well as the few 
notorious failures) of Landau Fermi liquid theory 3 make the subject an important area of modern 
condensed matter physics but one which we will not have time to explore. Instead, motivated 
in part by the success of Fermi liquid theory, we will proceed to explore the quantum partition 
function of the weakly interacting electron gas, r s <l. 



To prepare for a discussion of the field integral, we must first recast the Hamiltonian 
in second quantised form as 

H = J d Z rcl^-c a + \ J d'rd s r'ci(r)cl(r')j-^c a/ (r')c a (r), 

where the sum over repeated spin indicies a is here implied. When cast as functional field 
integral, the corresponding quantum partition function takes the form Z = tr(e _,s (' ff_/lJV )) = 
/ D(vjj,tl))e~ s \ where 




Employing the "four-vector" shorthand x = (r, r) and q = (u m ,q), where u m denotes 
a bosonic Matsubara frequency (exercise: think why), we may define the "density field" 
P q = -jpj dxe iqx Vv^VvOe) = -^Y.pi'pai'p+qa (with J dx = f^dr J d 3 r and q ■ x = 
UmT — q • r). Then, expressed in the Fourier basis, the action takes the form (exercise) 



2 As a simple non-interacting example, consider the adiabatic evolution of the bound states of a quan- 
tum particle as the confining potential is changed from a box to a harmonic potential well (sec Fig. 5.1). 
While the wavcfunctions and energies evolve, the topological characteristics of the wavef unctions, i.e. 
the number of the nodes, and therefore the assignment of the corresponding quantum numbers remains 
unchanged. 

3 L. D. Landau, Sov. Phys. JETP 3, 920 (1956); ibid. 5, 101 (1957). 
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where V(q) = / d 3 re~ l(l ' r V(r) = 47re 2 /q 2 . 4 Here, the prime on the q summation de- 
notes the exclusion of the q = contribution and reflects the presence of a neutralising 
background charge. 

Now, being quartic in the fields ip a , the Coulomb interaction prevents an explicit 
computation of the functional ^-integral. However, it is actually a straightforward matter 
to reduce, or 'decouple' the interaction operator bringing it to a form quadratic in the 
fields ip. Let us multiply the functional integral by the 'fat unity' 



1 = M / D(j) exp 



where <fi represents a complex bosonic field variable, and a normalization constant has 
been absorbed in the definition of the functional measuere D<p. Employing the variable 
shift e<p q I— > e<p q + jjj2V(q)p q , one obtains 



D(p exp 



L 3 /2 



+ 



1 



2L 3 



The rational behind this exercise can be seen in the last contribution to the exponent: This 
term is equivalent to the quartic interaction contribution to the fermionic path integral, 
albeit with opposite sign. Therefore, multiplication of Z by our unity leads to the field 
integral Z — J D<p f D(jp U) ip a )e~ s , where 



S 



-E 



-iui n H 

2m 



\x 5 pp/ + 



- : e(j)(p' — p) 



Ipp'a, (5.2) 



denotes the action, i.e. an expression that is free of quartic field interactions of ip a . Before 
proceeding, to acquire some intuition for the nature of the action, it is helpful to rewrite S 
in a real space/time representation. With 4> q = ^/=| / dx e~ tq ' x 4>(x) (four-vector notation 

as above), one may confirm that, 



dr / d r 



57T 



;(«90) 2 + ^ 



2m 



Physically, <fi couples to the electron degrees of freedom as a space/time dependent (imag- 
inary) scaler potential, while the first term reflects the Lagrangian energy density asso- 
ciated with the electric component of the electromagnetic (alias the photon) field. Said 
differently, the field represents the gauge particle — in this case, the photon — that 
mediates the Coulomb interaction between electrons. Before proceeding, let us now step 
back and discuss the general philosophy of the manipulations that led from the original 
partition function to the two-field representation (5.2). 

> INFO. The sequence of manipulations developed above, i.e. the 'decoupling' of a quartic 
interaction through an auxiliary field, is known more generally as a Hubbard-Stratonovich 



4 By representing Poisson's equation for a point charge in the Fourier space, one may confirm that 
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transformation. The essence of the transformation is a straightforward manipulation of a 
Gaussian integral. To make this point more transparent, let us reformulate the Hubbard- 
Stratonovich transformation in a notation that is not burdened by the presence of model-specific 
constants. Consider an interaction operator of the form Si nt = V a ^s ^aipp'^'yips (summation con- 
vention implied), where ip and ip may be either bosonic or fermionic field variables, the indices 
a, (3, . . . refer to an unspecified set of quantum numbers, Matsubara frequencies, etc., and V a p^s 
is an interaction matrix element. Now, let us introduce composite operators p a p = vpaipp to 
rewrite the interaction as S- mt = V a ^ys PapP-yS- The notation can be compactified still further 
by introducing composite indices m = (a/3), n = ("fd), whereupon the action Si nt = p m V mn p n 
acquires the structure of a generalized bilinear form. To reduce the action to a form quadratic 
in the ips one may simply multiply the exponentiated action by unity, viz. 

.-^./^.-^♦..-ww., 

" * ' 

1 

where <f> is bosonic. (Notice that here represents the matrix elements of the inverse and 
not the inverse {V mn )~ l of individual matrix elements.) Finally, applying the variable change 
<p m — ► cf) m + 2i(Vp) m where the notation (Vp) is shorthand for V mn p n , one obtains 



exp [-pmV mn p n ] = J " D(j) exp 



I.e. the term quadratic in p is cancelled. 5 This completes the formulation of the Hubbard- 
Stratonovich transformation. The interaction operator has been traded for (a) an integration 
over an auxiliary field coupled (b) to a V'-bilinear (the operator 4> m Pm)- 

> In essence, the Hubbard-Stratonovich transformation is tantamount to Gaussian integral 
identity (3.13) but read in reverse. An exponentiated square is removed in exchange for a 
linear coupling. (In (3.13) we showed how terms linear in the integration variable can be 
removed.) 

> To make the skeleton outlined above a well defined prescription, one has to be more 
specific about the meaning of the Gaussian integration over the kernel m V^0„, i.e. the 
integration variables can be real or complex, and V must be a positive matrix (which is 
usually the physical situation). 

> There is some freedom as to the choice of the integration variable. For example, the factor 
of 1/4 in front of the Gaussan weight (pmV^'Pn was introduced for mere convenience (viz. 
to generate a coupling (f) m p m free of numerical factors). If one does not like to invert the 
matrix kernel V mn , one can scale (ft m — > (V(j)) m , whereupon the key formula reads 



D(f>e 



> Exercise. Show that the passage from the Lagrangian to the Hamiltonian formulation of 
the Feynman path integral can be interpreted as a Hubbard-Stratonovich transformation. 



5 Here we have assumed that the matrix V is symmetric. If it is not, we can apply the relation 
PmYmnPn = p T Vp — | [/3 T (F + ^ T )/5] to symmetrize the interaction. 
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Figure 5.2: On the different channels of decoupling an interaction by Hubbard-Stratonovich 
transformation. Left: decoupling in the 'density' channel; middle: decoupling in the 'pairing' or 
'Cooper' channel; and right: decoupling in the 'exchange' channel. 

As defined, the Hubbard-Stratonovich transformation is exact. However, to make it a mean- 
ingful operation, it must be motivated by some physical considerations. In our discussion above, 
we split up the interaction by choosing p a p as a composite operator. However, there is clearly 
some arbitrariness with this choice. Why not, for example, pair the fermion-bilinears according 
to (t/j a ips){'4>'y'4>5), ° r otherwise? The three inequivalent choices of pairing up operators are shown 
in Fig. 5.2 where, as usual, the wavy line with attached field vertices represents the interaction, 
and the dashed ovals indicate how the field operators are paired. 

The version of the transformation discussed above corresponds to the left diagram of the 
figure. That type of pairing is sometimes referred to as decoupling in the direct channel. The 
denotation becomes more transparent if we consider the example of the spinfull electron-electron 
interaction, 

S'int = ^ J ' dr J d 3 rdV^ (T (r,r)^ (T /(r / ,r)F(r - r')Vv(r', r)^ CT (r, r), 

i.e. here a = f3 = (r, r, a), 7 = 6 = (r',r, a'), and V a pjg = V(r — r'). The 'direct' decoupling 
proceeds via the most obvious choice, i.e. the density operator p(r,r) = tp a (r,T)ip a (r,r). One 
speeks about decoupling in a 'channel' because, as will be elucidated below, the propagator of 
the decoupling field can be interpreted in terms of two Green function lines tied together by 
multiple interactions, a sequential object reminiscent of a 'channel'. 

However, more important than the terminology is the fact that there are other choices for 
p. Decoupling in the exchange channel is generated by the choice p ai ~ tpa^S where, in 
the context of the Coulomb interaction, the reversed pairing of field operators is reminiscent of 
the exchange contraction generating Fock-type contributions. Finally, one may decouple in the 
Cooper channel, p = V> a ^ 7 , p/3 7 = pj^g. Here, the pairing field is conjugate to two creation 
operators. Below we will see that this type of decoupling is tailored to problems involving 
superconductivity. 

The remarks above may convey the impression of a certain arbitrariness inherent in the 
Hubbard-Stratonovich scheme. Indeed, the 'correct' choice of decoupling can only be motivated 
by physical reasoning, not by plain mathematics. Put differently, the transformation as such is 
exact, no matter what channel we choose. However, later, we will want to derive an effective 
low energy theory based on the decoupling field. In cases where one has accidentally decoupled 
in an 'unphysical' channel, it will be difficult, if not impossible to distill a meaningful low 
energy theory for the field <j> conjugate to p. Although the initial model still contains the full 
microscopic information (by virtue of the exactness of the transformation) it is not amenable to 
further approximation schemes. 

In fact, one is frequently confronted with situations where more than one Hubbard-Stratonovich 
field is needed to capture the full physics of the problem. To appreciate this point, consider the 
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Coulomb interaction in momentum space. 

(5.3) 

Pl,...,P4 

In principle, we can decouple this interaction in any one of the three channels discussed above. 
However, 'interesting' physics is usually generated by processes where one of the three unbound 
momenta entering the interaction vertex is small. Only these interaction processes have a chance 
to accumulate an overall collective excitation of low energy (cf . many of the examples to follow) . 
It may be instructive to imagine the situation geometrically: In the three dimensional cartesian 
space of free momentum coordinates (pi,P2,P3) entering the vertex, there are three thin layers 
where one of the momenta is small, (q,P2,P3), (pi, Q,P3), (pi,P2i <l), \q\ "C \pi\- (Why not make all 
momenta small? Because that would be in conflict with the condition that the Green functions 
connecting to the vertex be close to the Fermi surface.) One will thus often choose to break down 
the full momentum summation to a restricted summation over the small-momentum sublayers: 

SintbP,^} - (^P^P+9^(q)^WVVp'-<7 - $<T P ^a>p+qV(p' ~ p)4>a'p>+q^ap> ~ 

p,p',q 

Vv P Vv_p +f ,y(p' - p)^ CT? /Vv-p'+g) • 

Now, each of these three contributions has its own predestined choice of a slow decoupling field. 
The first term should be decoupled in the direct channel f>d,q ~ Sp^opYVp+gj the second in 
the exchange channel p^aa'q ~ "ipap^^'p+q, and the third in the Cooper channel p c ,aa'q ~ 

VVpVv-p+q- One thus winds up with an effective theory that contains three independent 
slow Hubbard-Stratonovich fields. (Notice that the decoupling fields in the exchange and in the 
Cooper channel explicitly carry a spin-structure.) 

After this digression on the principles of the Hubbard-Stratonovich transformation, let us 
now return to the discussion of the electron gas. 



At the expense of introducing a second field, the Hubbard-Stratonovich transformation 
provides an action quadratic in the fermion fields. In this case, the fermion integration 
can be undertaken exactly. Making use of the Gaussian functional integral (4.19), one 
obtains 



= J D0 e -^ E «^ q2 ^ det 



\ P 2 

o T + fi + iecj) 

2m 



where the factor of two accounts for the spin degeneracy. 

The standard procedure to deal with the determinants generated at intermediate stages 
of the manipulation of a field integral is to simply re-exponentiate them. This is achieved 
by virtue of the identity, 



In det A = tr In A 



(5.4) 



valid for arbitrary (non-singular) operators A 6 Thus, the quantum partition function 



6 Eq. (5.4) is readily established by switching to an eigenbasis whereupon one obtains In det A = 
J2 a l n e a = tr In A, where e a are the eigenvalues of A and we have used the fact that the eigenvalues of 
In A are In e a . 
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takes the form Z — J D(pe where 

5[0] = ^^M 2 0- 9 -2tr In 



8tt 
q 



\ P 2 

o T + jj, + ie<j) 

2m 



(5.5) 



This is as far as purely formal exact manipulations can carry us. We have managed to 
trade the integration over the interacting Grassmann field ip a for an integration over an 
auxiliary field 0; a field that we believe encapsulates the relevant degrees of freedom of 
the model. This completes steps 1, 2, and 3 of the general programme outlined above. 

Ordinarily, the next step in the programme is to subject the action to a stationary 
phase analysis, i.e. to seek solutions of the set of saddle-point equations such that 

V 8 = (q#0, W ): ^1 = 0. 

b(p q 

Such a solution 0(x, t) <-> <p q is commonly referred to as a mean-field. This terminology 
can be understood by inspection of the argument of the 'tr In' above. The structure 
p 2 /2m — n + ie(f), where is a fixed configuration (to be determined by solving the 
saddle-point equations), resembles the Hamiltonian operator of particles subject to some 
background potential, or 'mean' field. The notation on the left hand side of the saddle- 
point equations indicates that our original interaction V(q) and, therefore, the decoupling 
field do not possess a zero momentum mode (a consequence of charge neutrality). 

However, in the present case, since the interaction is considered weak, we may an- 
ticipate that the solution to the saddle-point variation is the trivial one, viz. = 
- an assumption that we may check self-consistently. In this case, step 4 of the gen- 
eral programme may be considered as achieved and we may turn to explore fluctuations 
around = 0. Since the mean- field solution vanishes, it makes no sense to introduce new 
notation, i.e. we will denote the fluctuations again by the symbol 0. 

As regards the first term in the action (5.5), it has already a quadratic form. The 
logarithmic contribution can be expanded as if we were dealing with a function (again, a 
consequence of the trace), i.e. setting 

&- 1 = G 1 + i&j>, 

Gq 1 = d T + ^ — fM is the momentum and frequency diagonal operator whose matrix 
elements give the free Green function of the electron gas, we may express 

tr In G- 1 = tr ln(Go 1 + ie<j>) = tr In G 1 + tr ln(l + ieG (f)) 



2 



= tr In Gq 1 + ietr (G o 0) + -tr (G o 0G o 0) + . . . . 

Being 0-independent, the first term generates an overall constant multiplying the func- 
tional integral, viz. a constant that must describe the non-interacting content of the 
theory. Indeed, one may note that e 2tvXn6 ° l = e~ 2trlnda = (detC^ 1 ) 2 = Z is just the 
partition function of the non-interacting electron gas. Linear in 0, the second term of the 
expansion must, by virtue of the mean-field analysis, vanish. (Afterall, we are expanding 
around an extremum! To this end, one may note that tr (G o 0) = J2 q Go(q)4> g =o = 0.) 
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The third term is the interesting one. Remembering that couples to the theory as a 
potential, this term describes how potential flucutations are affected by the presence of 
the electron gas, i.e. it must encode the screening of the electromagnetic field by the 
electron degrees of freedom. 

To resolve this connection, let us make the momentum dependence of the second-order 
term explicit (exercise 7 ): 



where, setting £ p = ^ - /i, 



n(g) = ^G (p)G (p + g) = ^^ 



1 



1 



-iu n + £ p -iu n - iu m + £ 



p+q 



Collecting together the bare photon action with this expansion, to leading order in e, 
the quantum partition function takes the form 



-si*] 



Z = Z D<pe 



where the effective action for the electromagnetic field is given by 




This result has a clear physical interpretation: the interaction of the electron gas with 
the electromagnetic field induces a modified or screened Coulomb interaction (see 
Fig. 5.3), 



1 4tt 
e(uj m , q) q 2 ' 



e(w n ,q) = 1 —U(u m ,q) 



where e(u m ,q) is the energy and momentum dependent effective dielectric function. 
This perturbative result, which is known in the literature as the Random Phase Ap- 
proximation (RPA), amounts to treating the long-range part of the Coulomb inter- 
action as an "external" polarisation field, and the correction to the dielectric function, 
(47re 2 /q 2 )II(u; m , q), is known as the screened polarisability. To explore the dielectric 
properties of the interacting electron gas it is necessary to understand the frequency and 
momentum dependence of the density-density response function (5.6). To do so we must 
first learn how to perform Matsubara frequency summations. 



> INFO. Frequently, in working with imaginary time field integrals one often needs to 
perform Matsubara frequency summations. At low temperatures it can be argued that 



7 Hint: Consider the strategic incorporation of the resolution of identity 1 = J2 q underneath the 

trace. 
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k+q, co n +co n 
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Figure 5.3: The modified screened Coulomb interaction, D(u m ,q) can be viewed as the 
summation of an infinite 'diagrammatic' series expansion in the interaction: The bare 
Coulomb interaction vertex is 'dressed' by repeated particle-hole excitations of the electron 
gas. The corresponding summation of the geometric series (known as a Dyson equation) 
is shown schematically and can be compared to Eq. (5.6). 



the spacing 2tt/P between neighbouring Matsubara frequencies scales to zero which legitimates 
changing the summations into integrals. However, as we have seen previously, there exists an 
alternative way of evaluating the sums which not only keeps the finite temperature content of 
the theory but is also more efficient even in the limit of zero temperature. 

Referring to our ealier discussion for details, the basic idea behind the standard scheme of 
evaluating frequency summations is to introduce a complex auxiliary function g{z) that has 
simple poles at z = iu n . The sum f(iu> n ) one wishes to compute then emerges as the sum of 
residues obtained by integrating the product gf along a suitably chosen path C in the complex 
plane. The choice of both the function g and the integration contour C depend on the structure 
of the sum in which one is interested (convergence and analyticity properties of /, etc.). In the 
present case, the screened polarisability (5.6), 



9(z) 



exp(pz) + 1 



f(z) = 



z — li + e z — li + e' ' 



and C is a circular contour of infinite radius in the complex plane (see Fig. 5.4). We are therefore 
led to consider the integral 



1= f dz 
Jc 



1 



1 



q exp(f3z) + 1 z — li + e z — ll + e' 



(5.6) 



Two important observations can be made without explicit computation: (i) The integral exists 
(the integrand decays sufficiently fast in all directions for \z\ — > oo), and (ii) 1 = (the reason 
being that for \z\ — ► oo, the product fg < z~ J , where 7 > 1). Thus the integral of fg over a 
circle (radius ~ z) scales to zero as the radius — > 00. On the other hand, the integral along C 
gives the sum over the residues of all enclosed poles. The function fg has poles at the Matsubara 
frequencies iuj n and poles on the real axis at z = fx — e and z = ll — e' . Hence, 

= J^dz(fg)(z) = 2iri ^s (fg)(iu n ) + Res (fg)(» - e) + Res {fg)(» - e'^J . 

At the Matsubara frequencies, / is analytic and g has, by construction, unit residue. Thus 

Res (fg)(iu n ) = -/fan)- 
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Poles at Klatsubara 
frequencies 




Figure 5.4: Complex integration contour employed in calculating the sum (5.6). 

At z = \x — e, / has a simple pole, whereas g is analytic, i.e. 

Res (fg)^- e ) = 9 j^-. 
Combined with the pole at z = \x — e', and making use of the identity e 1 ^" = 1, we obtain 

TJ (l , „\- 2 ^ TO f(&) -M^k+q) 

LL{uj m ,q — —tt > — — 



(5.7) 



where we have made use of the fact that, on the real axis, the auxiliary function g is proportional 
to the Fermi-Dirac distribution function, 



^F(Ck) 



1 

e^(/i-e k ) + i ' 



(5.8) 



To analyse the screened Coulomb interaction D(u n ,q) we consider Eq. (5.7) for the 
density-density response function, and divide our consideration into two limits: 

> Static Limit (\u n \ <C k F \q\/m) 

For frequencies small as compared to the momentum transfer (and temperatures 
T <C fi), the response function converges to the static limit 

n(0,q) ^ q-d k n F (a) [d 3 k [°° 

where here we have made use of the fact that £k depends only on |k|, and we have 
deployed the continuum limit, i— > J = J °° dev{e) with z/(e) = l/|Vk£k| = 

2mA;/ (27r) 2 the density of states. From this result, one obtains the screened Coulomb 
interaction 

D(q) ' 



47T 
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When transformed back to real space, one obtains the effective static screened 
Coulomb interaction 



D(r) = -^ e -l r l/ ATF 



where Atf = (87re 2 z/(/i)) 1//2 defines the Thomas-Fermi screening length. Physi- 
cally, the bare Coulomb interaction is screened by the collective fluctuations of the 
electron gas. 

> High Frequency Limit {\uj n \ ^> kp\c{\/m) 

By contrast, in the high frequency limit, the density-density response function takes 
the form 

IL(u m , q) ~ -— \ — - ~ - / — — — — 1 + q • dkn F (&) 

L A ^ iuj m — q ■ <9 k 4 k J (2vr) 3 iu m \ imuJ m J 

by parts f d 3 k 2q 2 1 4 , 3 V N 2(f TIC? 



f d 6 k 2q 2 14, 

/ 7^3 Y Uf ^ k = — To T3 <j 

J (27r)- i mojf n (27r) d 3 



F muJm 2L 3 mu^ mu^ 

where n = N/ L d denotes the total number density. Applying the analytic continu- 
ation to real frequencies, iu) n — > u + iO, we obtain 



4:7ie 2 

lim D(u, q) = — — 

mw/fe_F|q|— +oo q 



0^ 

a; 2 



where w p = 47re 2 n/m represents the Plasma frequency. At high frequencies 
the dielectric response of the system is sensitive to the plasma oscillations of the 
background electron charge. In particular, for uo = u p , the collective excitations 
become undamped. 

> Ground State Energy 

Finally, from the partition function, it is possible to obtain an estimate of the ground 
state energy of the interacting electron gas. 

lim Z ~ e~ pEgs -. 

In the RPA approximation, performing the functional integral over the Gaussian 
action in 0, we obtain 

£ g .s. =£ g . B .(e = 0)-^5^]nD(w n ,q). 

where E g s _ = —1ylZq/(3 = 3n/i/5 is the ground state energy of the free electron gas. 
This was the formula derived by Gell-Mann and Bruckner from which, after some 
extensive algebra (exercise!), one obtains the high density (r s <C 1) expansion 

E g . s . = n ( ^ - + 0.622 In r s - 0.142 J Ryd. 

\ r s r s J 
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This concludes our preliminary analysis of the screening properties of the weakly 
interacting electron gas. By employing the coherent state path integral, we were 
(implicitly) able to establish the stability of the non-interacting ground state, deter- 
mine the modified screened Coulomb interaction, and identify plasma oscillations. 
In the following section, we will use the path integral to study the weakly interacting 
Bose gas and the phenomenon of superfluidity. 

5.3 Bose-Einstein Condensation and Superfluidity 

Previously, we have considered the influence of weak Coulomb interaction on the proper- 
ties of the electron gas. In the following, our goal will be to consider the phases realised 
by a weakly interacting Bose gas. To this end, let us introduce the quantum partition 
function Z = J D(rj>, ip) e~ s ^\ where 

<Kr, r)(d T + H - vMr, r) + |$(r, r)^(r, r)) 2 ] . (5.9) 

Here ip represents a complex field subject to the periodic boundary condition ip(r,(3) = 
ip(r, 0). The functional integral Z describes the physics of a system of bosonic particles in 
^-dimensions subject to a repulsive contact interaction of strength g > 0. For the moment 
the specific structure of the one-body operator H need not be specified. The most 
remarkable phenomena displayed by systems of this type are Bose-Einstein condensation 
and superfluidity. However, contrary to a widespread belief, these two effects do not 
mutually depend on each other: Superfluidity can arise without condensation and vice 
versa. We begin our discussion with the more elementary of the two phenomena. 

5.3.1 Bose-Einstein Condensation 

As may be recalled from elementary statistical mechanics, at sufficiently low temperatures, 
the ground state of a bosonic system can involve the condensation of a macroscopic 
fraction of particles into a single state. This phenomenon, predicted in a celebrated work 
by Einstein 8 is known as Bose-Einstein condensation. To see how this phenomenon is 
born out of the functional integral formalism, let us temporarily switch off the interaction 
and turn to the basis in which the one-particle Hamiltonian is diagonal. Expressed in the 
frequency representation, the partition function of the non-interacting system is given by, 

- ^2 ^ an (~ iUJ ™ + e « — Z 1 ) ^an ■ 

an 



Albert Einstein 1879-1955: 1921 Nobel Laureate in Physics "for his ser- 
vices to theoretical physics, and especially for his discovery of the law of 
the photoelectric effect". His work on the low temperature behaviour of 
the bosonic quantum gas is published in A. Einstein, Quantentheorie des 
einatomigen idealen Gases, Sitzungsber. Preuss. Akad. Wiss. 1925, 14 
8 (1925). 
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Without loss of generality, we may assume that the eigenvalues e a > are positive with 
a ground state €q = 0. 9 (In contrast to the fermionic systems discussed above, we should 
not have in mind low energy excitations superimposed on high energy microscopic degrees 
of freedom. Here, everything will take place in the vicinity of the ground state of the 
microscopic single-particle Hamiltonian.) Furthermore, we note that, to ensure stability, 
the chemical potential determining the number of particles in the system must be negative 
for, otherwise, the Gaussain weight corresponding to the low-lying states e a < — [i would 
change sign resulting in an ill-defined theory. 

From our discussion of section 4.2.1 we recall that the number of particles in the 
system is set by the relation (k B = 1) 



where, as usual, 713(e) = (e^ e ~^ — 1) _1 denotes the Bose distribution. For a given 
number of particles, this equation determines the temperature dependence of the chemical 
potential, (jl(T). As the temperature is reduced, the distribution function controlling 
the population of individual states decreases. Since the number of particles must be 
kept constant, this scaling must be counter-balanced by a corresponding increase in the 
chemical potential. 




Figure 5.5: Schematic plot showing the variation of the chemical potential as a function of 
temperature. Note that, for T < T c , the chemical potential remains pinned at zero. 

Below a certain critical temperature T c , even the maximum value of the chemical 
potential, [i = 0, will not suffice to keep the distribution function 7i B (e a ^ ) large enough to 

accomodate all particles in the states of non-vanishing energy, viz. ^2 a>0 TiB(e a )| M=0 =° 
Ni < N. I.e. below the critical temperature, the chemical potential stays constant at 
H = (see the figure). As a result, a macroscopic number of particles N — Ni must 
accumulate in the single-particle ground state: Bose— Einstein condensation. 

> EXERCISE. For a three-dimensional free particle spectrum, e& = h 2 k 2 /2m, show that the 
critical temperature is set by T c = ^3-, where a = p -1 / 3 is the average interparticle spacing, 
and Co is a constant of order unity. Show that for temperatures T < T c , the density of particles 
in the condensate (k = 0) is given by po(T) = p[l — (yp-) 3 ' 2 ]- 



> INFO. Since its prediction in the early 20s, the phenomenon of Bose-Einstein condensation 
has been a standard component of undergraduate texts. However, it took some seven decades 

9 The chemical potential /i can always be adjusted so as to meet this condition. 
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before the condensation of bosonic particles was directly observed in experiment. The reason 
for this delay is that the critical condensation temperature of particles that are comfortably 
accessible to experiment — atoms — is absurdly low. 

In 1995 the groups of Cornell and Wieman at Colorado University and, soon after, Ketterle 
at MIT succeeded in cooling a system of rubidium atoms down to temperatures of 20 billionths(l) 
of a Kelvin. 11 To reach these temperatures, a gas of rubidium atoms was caught in a magnetic 
trap, i.e. a configuration of magnetic field gradients that couple to the magnetic moments of 
the atoms so as to keep the system spatially localized (see the schematic). 




Figure 5.6: Spectroscopic images of a gas of atoms at 400nK (left), 200nK (middle), and 50nK 
(right). The peak in the density distribution signals the onset of condensation. Courtesy of 
JILA institute, University of Colorado 

The gas of atoms was then brought to a temperature of O(10~ 5 )K — still much too hot 
to condense — by 'laser cooling'; crudely speaking, a technique where atoms, subjected to a 
suitably adjusted ray of monochromatic light, may transmit more of their kinetic energy to the 
photons than they get back. To lower the temperature still further, the principle of 'evaporative 
cooling' was applied: By lowering the potential well of the trap, a fraction of the atoms, namely 
those with large kinetic energy, is allowed to escape. The remaining atoms have a low kinetic 
energy and, therefore, a low temperature. What sounds like a simple recipe actually represents 
a most delicate experimental procedure. (For example, if the trap potential is lowered too 
strongly, all atoms escape and there is nothing left to condense. If, on the other hand, trapping 
is too strong, the atoms remain too hot, etc.) However, after more than a decade of intensive 
experimental preparation, the required temperatures have been reached. 

Spectroscopic images of the Bose-Einstein condensation process are shown in the figure 
(courtesy of the JILA institute, University of Colorado) for three values of temperature (400 
nK, 200 nK, and 50 nK from left to right). The peak in the density distribution signals the onset 
of condensation. On lowering the temperature, one may observe the transition to a condensed 
phase by monitoring the formation of a peak in the density distribution. The preparation of 

10 Here, by 'direct' we refer to the controlled preparation of a state of condensed massive bosonic 
particles. There are numerous 'indirect' manifestations of condensed states, e.g. the anomalous properties 
of the Helium liquids at low temperatures, or of Cooper-pair condensates in superconductors. 

n M. H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, and E. A. Cornell, Observation of 
Bose-Einstein Condensation in a Dilute Atomic Vapor, Science 269, 198 (1995). 
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a Bose condensed state of matter was recognized with the award of the 2001 Nobel prize in 
physics. Since 1995, research on atomic condensates has blossomed into a broad arena of re- 
search. Already, it is possible to prepare complex states of Bose condensed matter such as atomic 
vortices in rotating Bose-Einstein condensates, condensates in different dimensionalities, or even 
an artificial crystalline state of matter. A detailed discussion of these interesting developments 
is beyond the scope of the present text. 



With this background, let us now try to understand how the phenomenon of Bose- 
Einstein condensation can be implemented into the functional integral representation. 
Evidently, the characteristics of the condensate will be described by the zero field compo- 
nent ipoij). The problem with this zero mode is that, below the condensation transition, 
its action apears to be unbound: both the chemical potential and the eigenvalue are zero. 
This means that the action of the zero Matsubara component ^0,0 vanishes. We will deal 
with this difficulty in a pragmatic way. That is, we will not treat ipo(r) as an integration 
variable but rather as a time-independent Lagrange multiplier to be used to fix the num- 
ber of particles below the transition. More precisely, we introduce a reduced action of the 
form 

So[$0, fpo] = -$o/3^0 + ^2 ^ an ( iUJn + e a- ^an , 

where we did not yet set /i — (since we still need \i as a differentiation variable). To 
understand the rational behind this simplification one may note that 

N = -d M F | M=0 - = T^lnZ | M=0 - = ^0 + T ^ . \ = ^0 + N, (5.10) 

determines the number of particles. According to this expression, ipoipo = No sets the 
number of particles in the condensate. Now, what enables us to regard ipo as a time- 
independent field? Remembering the construction of the path integral, we note that 
the introduction of time-dependent fields, or 'time slicing' was necessitated by the fact 
that the operators appearing in the Hamiltonian of a quantum theory do not, in gen- 
eral, commute. (Otherwise we could have decoupled the expression tr [e~^ H ~ lxN ^ a '°^) ~ 
J d(*j),ip)e~ f3( - H ~^ N ^^'^ in terms of a single coherent state resolution, i.e. a 'static' con- 
figuration). Reading this observation in reverse, we conclude that the dynamic content of 
the field integral represents the quantum character of a theory. (Alluding to this fact, the 
temporal fluctuations of field variables are often referred to as quantum fluctuations.) 
Conversely, a static approximation in a field integral ip(r) = ipo = const, amounts to 
replacing a quantum degree of freedom by its classical approximation. 
(In order to distinguish them from quantum, fluctuations in the 'classical' static sector of 
the theory are called thermal fluctuations.) To justify the approximation of a <-> t/j 
by a classical object, notice that, upon condensation, Nq = (ajao) will assume 'macro- 
scopically large' values. On the other hand, the commutator [a > a o] = 1> continues to be 
of 0(1). It thus seems to be legitimate to neglect all commutators of the zero operator 
cio in comparison with its expectation value — a classical approximation. 12 

12 Notice the similarity of that reasoning to the arguments employed in connection with the semi- 
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Now, we are still left with the problem that the ^"integration appears to be un- 
defined. The way out is to remember that the partition function should extend over those 
states that contain an (average) number of iV particles. That is, Eq. (5.10) has to be 
interpreted as a relation that fixes the modulus ipoipo s° as to adjust the appropriate value 
of N. (For a more rigorous discussion of the choice of the thermodynamic variables in the 
present context, we again refer to Ref.[l].) 

5.3.2 The Weakly Interacting Bose Gas 

Now, with this background, let us restore the interaction focussing on a small but finite 
coupling constant g. To keep the discussion concrete, we specialize to the case of a free 
single-particle system, H = p 2 /2m. (Notice that the ground state wavefunction of this 
system describes a spatially constant zero momentum state.) By adiabatic continuity we 
expect that much of the picture developed above will survive generalization to non-zero 
interaction strengths. In particular, the ground state, which in the case under consider- 
ation corresponds to a temporally and spatially constant mode ipo, will continue to be 
macroscopically occupied. Under these circumstances, the dominant contribution to the 
action will again come from the classical ip sector: 

TSi^M = -MxA) + ^(^o) 2 • (5.11) 

Crucially, the stability of the action is now guaranteed by the interaction vertex, no matter 
how small is g > (see the schematic plot of the action in the figure). Accordingly, 
we will no longer treat ipo as a fixed parameter but rather as an ordinary intergration 
variable. Integration over all field components will produce a partition function Z{p) 
that depends parametrically on the chemical potential. As usual in statistial physics, the 
latter can then be employed to fix the particle number. (Notice that, vis-a-vis aspects 
of thermodynamics, the interacting system appears to behave more 'naturally' than its 
ideal, non-interacting approximation. This reflects a general feature of bosonic systems; 
interactions 'regularize' a number of pathological features of the ideal gas.) 

Returning to the -?/>o~integration, we observe that, for low enough temperatures, the 
problem is an ideal candidate for saddle-point analysis. Variation of the action with 
respect to ip Q obtains 




This equation is solved by any constant complex field configuration ipo with modulus 
l^o I — a/ I^L d /g = 7. Inspite of its innocent appearance, this equation reveals much 
about the nature of the system: 

classical treatment of spin systems in the limit of large S (section 2.2.4). Unfortunately, the actual 
state of affairs with the classical treatment of the condensate is somewhat more complex than the simple 
argument above suggests. (For a good discussion, see Ref. [1].) However, the net result of a more thorough 
analysis, i.e. an integration over all dynamically fluctuating components ?/>o,n^o> shows that the treatment 
of ipo as classical represents a legitimate approximation. 
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Re $ 

Figure 5.7: The action S[ipo, tpo] shown as a function of real and imaginary part of the 
condensate field (part removed for clarity). The most important features of the action 
are (a) the existence of a degenerate minimum determined by the set of solutions of the 
equation d\^ \S = 0, and (b) the large amplitude asymptotics ~ /9fi'|V'o| 4 stabilizing the 
^o-integration. 

> For n < (i.e. above the condensation threshold of the non-interacting system), 
the equation exhibits only the trivial solution ip = 0. This means that no stable 
condensate amplitude exists. 

> Below the condensation threshold (i.e. for \x > 0), 13 the equation is solved by any 
configuration with |t/> | — 7 = a/ fJ,L d /g. (Notice that V'oV'o ^ L d , reflecting the 
macroscopic population of the ground state.) 

> The equation couples only to the modulus of ipo- he. the solution of the stationary 
phase equation is continuously degenerate: Each configuration ip = j exp(i(fi) , (ft G 
[0, 2ir] is a solution. 

For our present discussion, the last of the three aspects mentioned above is the most 
important. It raises the question as to which of the configurations ip = 7exp(?0) is the 
'right' one? 

Without loss of generality, we may chose tpo = 7 G R as a reference configuration for 
our theory. This choice amounts to selecting a particular minimum lying in the 'mexican 
hat' profile of the action shown above. However, it is clear that an expansion of the action 
around that minimum will be singular: Fluctuations tp — > ip + 5ip that do not leave the 
azimuthally symmetric well of degenerate minima do not change the action and, therefore, 
have vanishing expansion coefficients. As a result, in the present situation, we will not be 
able to implement a simple scheme viz. 'saddle-point plus quadratic fluctuations'. (There 
is nothing that constrains the deviations 5ip to be small.) The integral over fluctuations 
around the mean-field configuration has to be undertaken in a more careful way. 

> INFO. The mechanism encoutered here is one of spontaneous symmetry breaking. 

To understand the general principle, consider an action S[ip] with a global continuous symmetry 

13 Due to the stabilization of the zero mode integration by the interaction constant, /i < is no longer 
a strict condition. 
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under some transformation g (not to be confused with the afore mentioned coupling constant 
of the Bose gas): Specifically, the action remains invariant under a global transformation of the 
fields such that, Vi G M : tpi — > gipi, where M is known as the "base manifold" , i.e. S[ip] = S[gtp]. 
The transformation is 'continuous' in the sense that g takes values in some manifold, typically 
a group G. 

Examples: The action of a Heisenberg ferromagnet is invariant under rotation of all spins si- 
multaneously by the same amount, Sj — > gSi. In this case, g £ G = 0(3), the three-dimensional 
group of rotations. The action of the displacement fields u describing elastic deformations of a 
solid (phonons) is invariant under simultaneous translation of all displacements u» — ► Uj + a, 
i.e. the symmetry manifold is the the d-dimensional translation group G — M. d . In the example 
above, we encountered a U(l) symmetry under phase multiplication — ► e'^o- This phase 
freedom expresses the global gauge symmetry of quantum mechanics under transformation 
by a phase, a point we will discuss in more detail below. 

Now, given a theory with globally G invariant action, two scenarios are conceivable: Ei- 
ther the ground states share the invariance properties of the action or they do not. The two 
alternatives are illustrated in the figure for the example of the Bose system. For fi > 0, the 
action S^OjV'o] has a single ground state at V'o = 0. This state is trivially symmetric under 
the action of G = U(l). However, for negative fi, i.e. in the situation discussed above, there 
is an entire manifold of degenerate ground states, defined through the relation \i/jq\ = 7. These 
ground states transform into each other under the action of the gauge group. However, none of 
them is individually invariant. 

With the other examples mentioned above, the situation is similar. For symmetry groups 
more complex than the one-dimensional manifold U(l), the ground states will, in general, be 
invariant under transformation by the elements of a certain subgroup H C G (that includes the 
two extremes H = {1} and H = G). For example, below the transition temperature, the ground 
state of the Heisenberg magnet will be given by (domainwise) aligned configurations of spins. 
Assuming that the spins are oriented along the z-direction, the ground state is invariant under 
the abelian subgroup H C 0(3) containing all rotations around the z-axis. However, invariance 
under the full rotation group is manifestly broken. Solids represent states where the translation 
symmetry is fully broken, i.e. all atoms collectively occupy a fixed pattern of spatial positions 
in space, H = {1}, etc. 

Inspite of the undeniable existence of solids, magnets, and Bose condensates of definite phase, 
the notion of a ground state that does not share the full symmetry of the theory may appear 
paradoxical, or at least 'unnatural'. For example, even if any particular ground state of the 
'Mexican hat' potential shown in the figure above 'breaks' the rotational symmetry, shouldn't 
all these states enter the partition sum with equal statistical weight, such that the net outcome 
of the theory is again fully symmetric? 

To understand why symmetry breaking is a 'natural' and observable phenomenon, it is 
instructive to perform a gedanken experiment: To this end, consider the partition function of a 
classical 14 ferromagnet, 



where H is the rotationally invariant part of the energy functional and h represents a weak 
external field. (Alternatively, we can think of h as an internal field, caused by a slight struc- 
tural imperfection of the system.) In the limit of vanishing field strength, the theory becomes 

14 The same argument can be formulated for the quantum magnet. 
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manifestly symmetric. Symbolically, 

lim lim Z — ► rot. sym. , 

where the limit N — > oo serves as a mnemonic indicating that we consider systems of macroscopic 
size. However, keeping in mind the fact that the model ought to describe a physical magnetic 
system, the order of limits taken above appears questionable. Since the external perturbation 
couples to a macroscopic number of spins, a more natural description of an 'almost' symmetric 
situation would be 

lim lim Z — ► ? 

The point is that the two orders of limits lead to different results. In the latter case, for any h, 
the N — ► oo system is described by an explicitly symmetry broken action. No matter how small 
the magnetic field, the energetic cost to rotate N — ► oo spins against the field is too high, i.e. the 
ground state |S) below the transition temperature will be uniquely aligned, Sj || h. When we 
then send h — > in a subsequent step, that particular state will remain the observable reference 
state of the system. Although, formally, a spontaneous thermal fluctuation rotating all spins 
by the same amount |S) — > \gS) would not cost energy, that flucutation can be excluded by 
entropic reasoning. 15 (By analogy, one rarely observes kettles crashing into the kichen wall as a 
consequence of a concerted thermal fluctuation of the water molecules!) 

However, the appearance of non-trivial ground states is just one manifestation of sponta- 
neous symmetry breaking. Equally important, residual fluctuations around the ground state 
lead to the formation of soft modes (sometimes known as massless modes), i.e. field config- 
urations q whose action S[(f>] vanishes in the limit of long wavelengths q — ► 0. Specifically, the 
soft modes formed on top of a symmetry broken ground state are called Goldstone modes. 
As a rule, the presence of soft modes in a continuum theory has important phenomenological 
consequences. To understand this point, notice that the general structure of a soft mode action 
is given by 

S[4>] = J> q + c 2^ 2 ] </>- q + O(0 4 , q 3 ) , (5-12) 

q,i 

where c\ 2 are coefficients. The absence of a constant contribution to the action (i.e. a contri- 
bution that does not vanish in the limit q — ► 0) signals the existence of long-ranged power-law 
correlations in the system. As we will see shortly, the vanishing of the action in the long wave- 
length limit q — > further implies that the contribution of the soft modes dominates practically 
all observable properties of the system. 



5.3.3 Superfluidity 

As we have seen, the theory of the weakly interacting superfluid to be discussed below 
was originally conceived by Bogoliubov, then in the language of second quantisation. 16 In 
the following, we will reformulate the theory in the language of the field integral starting 

15 Note that this (overly) simple picture in fact breaks down in dimensions d < 2, cf. our discussion of 
the thermal fluctuations of the fcrromagnct in chapter 2. 

16 Bogoliubov, N. N., On the Theory of Superfluidity, J. Phys. (USSR) 11, 23 (1947), (reprinted in D. 
Pines, The Many-body Problem, W. A. Benjamin, New York, 1961). 
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with the action of the weakly interacting Bose gas (5.9). Focussing on temperatures 
below T c , (p > 0), let us expand the theory around the particular mean-field ground 
state 4>o = ipo = (pL d /g) 1 / 2 = 7. (Of course, any other state lying in the 'Mexican 
hat' minimum of the action would be just as good.) Notice that the quantum ground 
state corresponding to the configuration if> is unconventional in the sense that it cannot 
have a definite particle number. The reason is that, according to the correspondence 
if) <-> a between coherent states and operators, respectively, a non-vanishing functional 
expectation value of ipo is equivalent to a non-vanishing quantum expectation value (do). 
Assuming that, at low temperatures, the thermal average (. . .) will project onto the 
ground state \ fl), we conclude that (O|a |fi) 7^ 0, i.e. \fl) cannot be a state with a definite 
number of particles. 17 

The symmetry group U(l) acts on this state by multiplication, if> — > e^ipo and if>o — > 
e~ % ^i\)Q. Knowing that the action of a weakly modulated field 0(r, r) will be massless, let 
us introduce coordinates 

" Im if) \ S = extr. 

5p 

*- — *- 

Re ifj 



Figure 5.8: Schematic diagram showing the coordinates of the massive (5p) and massless (5(fi) 
fluctuations. 



^(r,r) = [p + Mr,r)] 1 /V^), 
^(r,r) = [p + 5p(r,r)] 1 /V^, 

where po = 7 2 = ^o^o is the condensate density. Evidently, the variable 5p parameterizes 
deviations of the field ip(r,r) from the extremum. These excursions are energetically 
costly, i.e. 5p will turn out to be a massive mode. Also notice that the transformation of 
coordinates — > {Sp, 0), viewed as a change of integration variables, has a Jacobian 
of unity. 

> INFO. As we are dealing with a (functional) integral, there is a lot of freedom as to the 
choice of integration parameters. (I.e. in contrast to the operator formulation, there is no a 
priori constraint for a transform to be 'canonical'.) However, physically meaningful changes of 
representation will usually be canonical transformations, in the sense that the corresponding 

17 However, as usual with grand canonical descriptions, in the thermodynamic limit, the relative uncer- 
tainty in the number of particles, ((N 2 ) — (N) 2 ) / (N) 2 will become vanishingly small. 
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transformations of operators would conserve the commutation relations. Indeed, one may con- 
firm that the operator transformation a(r) = p(r) 1 / 2 e % ^ r \ a\r) = e~ l ^ r ^ /5(r) 1//2 , fulfills this 
criterion (exercise). 



We next substitute the density-phase relation into the action and expand to second 
order around the reference mean-field. Ignoring gradients acting on the density field (in 
comparison with the 'potential' cost of these fluctuations), we obtain 



S[6p, 0] « J d d r J dr 



zpd T d> + —(d<f>) + — 



(5.13) 



The first term of the action has the canonical structure 'momentum x d T (coordinate)' 
indicative of a canonically conjugate pair. The second term measures the energy cost of 
spatially varying phase flucutations. Notice that fluctuations with 0(r, r) = const, do 
not incur an energy cost — is a Goldstone mode. Finally, the third term records the 
energy cost of massive fluctuations from the potential minimum. Eq. (5.13) represents the 
Hamiltonian version of the action, i.e. an action comprising coordinates <fi and momenta 
5 p. Gaussian integration over the field 5p leads us to the Lagrangian form of the action 
(exercise): 



S[(j>] ^^Jd d rjdr 



g m 



(5.14) 



Comparison with Eq. (1.2) identifies this action as the familiar the <i-dimensional oscil- 
lator. Drawing on the results of chapter 1 (see, e.g., Eq. (1.15)), we find that the en- 
ergy c^k carried by elementary excitations of the system scales linearly with momentum, 
u k = \k\p /mg. 

Let us now discuss the physical ramifications of these results. The actions (5.13) and 
(5.14) describe the phenomenon of superfluidity. To make the connection between the 
fundamental degree of freedom of a superfluid system, the supercurrent, and the phase 
field explicit, let us consider the quantum mechanical current operator 

j(r,r) = J- [(Vat(r,r))a(r,r) -at(r,r)Va(r,r)] " 

- ^ [W(r, r))V(r, r) - ^(r, r)V^(r, r)] » ^V0(r, r), (5.15) 

where the arrow indicates the functional integral correspondence of the operator descrip- 
tion and we have neglected all contributions arising from spatial fluctuations of the density 
profile. (Indeed, these — massive — fluctuations describe the 'normal' contribution to 
the current flow.) 



> INFO. Superfluidity is one of the most counterintuitive and fascinating phenomena 
displayed by condensed matter systems. Experimentally, the most straightforward access to 
superfluid states of matter is provided by the Helium liquids. Representative of many other 
effects displayed by superfluid states of Helium, we mention the capability of thin films to flow 
up the walls of a vessel (if the reward is that on the outer side of the container a low lying basin 
can be reached — the fountain experiment) or to effortlessly propagate through porous media 
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that no normal fluid may penetrate. 



The gradient of the phase variable is therefore a measure of the (super) current flow 
in the system. The behaviour of that degree of freedom can be understood by inspection 
of the stationary phase equations — alias, the Hamilton or Lagrange equations of motion 
- associated with the actions (5.13) or (5.14). Turning to the Hamiltonian formulation, 
one obtains (exercise) 

id T <p = —gSp, id T 5p = —d 2 (f) = V • j. 

m 

The second of these equations represents (the Euclidean time version) of a continuity 
equation. A current flow with non-vanishing divergence is accompanied by dynamical 
distortions in the density profile. The first equation tells us that the system adjusts to 
spatial fluctuations of the density by a dynamical phase fluctuation. The most remarkable 
feature of these equations is that they possess steady state solutions with non-vanishing 
current flow. Setting d T <f) = d T 5p = 0, we obtain the conditions 5p = and V • j = 0, i.e. 
below the condensation temperature, a configuration with a uniform density profile can 
support a steady state divergenceless (super) current. Notice that a 'mass term' in the 
(j) action would spoil this property, i.e. within our present approach, the phenomenon of 
supercurrent flow is intimately linked to the Goldstone mode character of the (f> field. 

> EXERCISE. Add a ficticious mass term to the faction (viz. 5C = mcf 2 ) and explore its 
consequences. How do the features discussed above present themselves in the Lagrange picture? 

It is very instructive to interpret the phenomenology of supercurrent flow from a dif- 
ferent, more microscopic perspective. Steady state current flow in normal environments 
is prevented by the mechanism of energy dissipation, i.e. particles constituting the 
current flow scatter off imperfections inside the system thereby converting part of their 
energy into the creation of elementary excitations. (Macroscopically, the conversion of 
kinetic energy into the creation of excitations manifests itself as heat production.) Ap- 
parently, this mechanism is inactivated in superfluid states of matter, i.e. the current flow 
is dissipationless. 

How can the dissipative loss of energy be avoided. Trivially, no energy can be ex- 
changed if there are no elementary excitations to create. In reality, this means that the 
excitations of the system are energetically high-lying such that the kinetic energy stored 
in the current-carrying particles is insufficient to create them. But this is not the situ- 
tation that we encounter in the superfluid! As we saw above, there is no energy gap 
separating the quasi-particle excitations of the system from the ground state. Rather, 
the dispersion u(k) vanishes linearly as k — > 0. However, there is an ingenuous argu- 
ment due to Landau showing that a linear excitation spectrum indeed suffices to stabilize 
dissipationless transport: 

D> INFO. Consider the flow of some fluid through a pipe (cf. Fig. 5.9 top left). To be concrete, 
let us assume that the flow occurs at a uniform velocity V. Taking the mass (of a certain portion 
of the fluid) to be M, the current carries a total kinetic energy E\ = MV 2 /2. Now, suppose we 
view the situation from the point of view of the fluid, i.e. we perform a Galileian transformation 
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into its own rest frame (see Fig. 5.9, top right). From the perspective of the fluid, the walls of 
the pipe appear as though they were moving with velocity —V. Now, suppose that frictional 
forces between fluid and the wall lead to the creation of an elementary excitation of momentum 
p and energy e(p), i.e. the fluid is no longer at rest but carries kinetic energy. After a Galileian 
transformation back to the laboratory frame, one finds that the energy of the fluid after the 
creation of the excitation is given by (excercise) 

MV 2 

E 2 = + p • V + e(p) . 

Now, since all of the energy needed to manufacture the excitation must have been provided by 





Figure 5.9: Top left: Flow of a fluid through a rough pipe. Top right: The same viewed from 
the rest frame of the fluid. Bottom left: Dissipative creation of a (quasi-particle) excitation. 
Bottom right: The same viewed from the laboratory frame. 

the liquid itself, energy conservation requires that E\ = E 2 , or — p • V = e(p). Since p • V > 
— |p||V|, this condition can only be met if |p||V| > e(p). While systems with a 'normal' gapless 
dispersion, e(p) ~ p 2 are compatible with this energy-balance relation (i.e. no matter how 
small |V|, quasi-particles of low momentum can always be excited), both gapped dispersions 
e(p) - — » const, and linear dispersions are incompatible if V becomes smaller than a certain 
critical velocity V*. Specifically for a linear dispersion e(p) = v\p\, the critical velocity is 
given by V* = v. For currents slower than that, the flow is necessarily dissipationless. 



Let us conclude our preliminary discussion of the weakly interacting Bose gas with a 
very important remark. Superficially, Eqs. (5.13) and (5.14) suggest that we have managed 
to describe the long-range behaviour of the condensed matter system in terms of a free 
Gaussian theory. However, one must recall that is a phase field, defined only modulo 
27T. (In Eqs. (5.13) and (5.14) this condition is understood implicitly. At this point, it 
is perhaps worth reiterating that when dealing with Goldstone modes it is important to 
keep the underlying geometry in mind and not too tightly focus on a specific coordinate 
representation.) The fact that is defined only up to integer multiples of 2ir manifests 
itself in the formation of the most interesting excitations of the superfluid; vortices, i.e. 
phase configurations <p(r, r) that change by a multiple of 2% as one moves around a certain 
reference coordinate, the vortex centre. Existing in parallel with harmonic phonon-like 
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excitations discussed above, these excitations lead to a wealth of observable phenomena. 
However, leaving such effects aside, let us turn to the discussion of another prominent 
superfluid, the condensate of Cooper pairs, more generally known as the superconductor. 



5.4 Superconductivity 

The electrical resistivity of many metals and alloys drops suddenly to zero when the 
specimen is cooled to a sufficiently low temperature. This phenomenon, which goes by 
the name superconductivity, was first observed by Kammerlingh Onnes in Leiden in 
1911, three years after he first liquefied Helium. 18 Perhaps more striking, a superconductor 
cooled below its transition temperature in a magnetic field expels all magnetic flux from 
its interior. This phenomenon of perfect diamagnetism is known as the Meissner effect 
and is characteristic of superconductivity. 




Figure 5.10: The levitation of a magnet above a (high temperature) superconductor due 
to the expulsion of magnetic flux. 



The superconducting state is an ordered state of the conduction electrons of the metal. 
The nature and origin of the ordering was explained by Bardeen, Cooper and Schrieffer. 19 
At low temperatures, the presence of an attractive pairwise interaction can induce an 
instability of the electron gas towards the formation of bound pairs of time-reversed states 

18 

Kammerlingh Onnes 1853-1926 
(left, photographed with van der 
Waals): 1913 Nobel Laureate 
in Physics for his investigations 
on the properties of matter at 
low temperatures which led, in- 
ter alia to the production of liq- 
uid helium. 

19 J. Bardeen, L. N. Cooper and J. R. Schrieffer, Phys. Rev. 106, 162 (1957); 108, 1175 (1957). 

John Bardeen 1908-1991 (left), Leon N. 
Cooper 1930- (centre), and J. Robert 
Schrieffer 1931- (right): 1972 Nobel Lau- 
reate in Physics for their jointly developed 
theory of superconductivity. (Bardeen 
was also recipient of the 1956 Nobel Lau- 
reate in Physics for his research on semi- 
conductors and discovery of the transistor 
effect. 
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k | and — k 1 in the vicinity of the Fermi surface. We have already seen (in section 2.2.6) 
how the exchange of lattice vibrations or phonons can induce an attractive interaction 
of electrons within the Debye frequency ujd of the Fermi surface. Being made up of two 
electrons, these composite objects, known as Cooper pairs behave as bosons. At low 
temperatures, these bosonic degrees of freedom form a condensate which is responsible 
for the remarkable properties of superconductors such as perfect diamagnetism. 

To explore the phenomenology of the Cooper instability of the electron gas, we will 
adopt a simplified model known as the pairing or reduced Hamiltonian 



Although, strictly speaking, a realistic model of attraction would involve a more compli- 
cated momentum-dependent interaction such as the one obtained from the consideration 
of the electron-phonon interaction in section 2.2.6, the simple pairing interaction cap- 
tures the essential physics. More importantly, to simplify our discussion, we will take the 
electrons to be otherwise non-interacting. In fact, the presence of a repulsive Coulomb 
interaction of the electrons plays a crucial role in the controlling properties of the super- 
conductor, a point to which we will return later. 



Figure 5.11: Schematic diagram showing the Fermi surface of the electron gas. The 
attractive interaction mediated by the exchange of phonons allows electrons within the 
Debye frequency of the Fermi surface to pair. 

Before turning to the field theoretic formulation, we will begin by investigating a 
mean-field theory of the BCS transition from the stand point of the second quantisation. 

5.4.1 Mean-Field Theory of Superconductivity 

In its present form, the Hamiltonian explicitly involves a two-body interaction of the 
electrons. As such, it seems infeasible to develop an exact many-body treatment of the 





Concepts in Theoretical Physics 



172 CHAPTER 5. BROKEN SYMMETRY AND COLLECTIVE PHENOMENA 



Hamiltonian. Instead we will seek an approximation which renders the Hamiltonian bi- 
linear in the electron operators and, therefore, tractable. As usual in physics, our method 
relies on the expected structure of the ground state wavefunction. In particular, antici- 
pating that electrons in time-reversed states pair, let us suppose that 

A = #^(g.s.|c_ ki c kT |g.s.), A* = J^(g.s.|cj cT c! k Jg.s.) 



acquires a non-zero expectation value in the ground state. Here A represents an order 
parameter becoming non-zero in the condensed phase and therefore signalling the tran- 
sition to the superconducting state. At first sight, a non-zero expectation value A looks 
strange: such a result would imply that the ground state wavefunction of the supercon- 
ducting condensate is not an eigenstate of particle number (while one can see that the 
Hamiltonian commutes with N). However, later, we will see that in the grand canoni- 
cal ensemble the ground state wavefunction is a superposition of states involving many 
particles but strongly peaked around the thermodynamic density N/L d . 
To develop the mean-field approximation, let us set 



small 



kJ. C kT 



A + UE c -k|C kT -A 



and keep only terms which depend up to quadratic order in the electron operators. Adding 
the chemical potential, the 'mean-field' Hamiltonian takes the form 



9 



known as the Bogoluibov or Gor'kov Hamiltonian. In this simplified form, it is 
interesting to note that the Hamiltonian does not now conserve particle number. Instead, 
pairs of particles are born and annihilated out of the vacuum. 

To bring the mean-field Hamiltonian to a diagonal form, it is convenient to recast it 
in a Nambu spinor representation defining 



t 



J 

-kT 



-kT 

t 

-kj. 



after which the Hamiltonian takes the form (exercise: recall the fermionic anticommuta- 
tion relations of the electron operators) 




Now, being bilinear in the electron operators, the mean-field Hamiltonian can be brought 
to a diagonal form by employing the unitary transformation 



a 



kT 



a 



-kT 



COS#k 

sin 6> k 



sin 6>k 
- cos 9^ 



^kT 

: -kT 



k' 
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(under which the anticommutation relations of the new electron operators a ka are main- 
tained: exercise). Note that the notation is purely symbolic: involves a superpo- 
sition of cL and c_kj- Choosing A to be real, 20 and setting tan(2#k) = — A/£k, i-e. 
cos(2# k ) = ^k/Ak, sin(2# k ) = — A/A k , where A k = (A 2 + £k) 1//2 > the transformed Hamilto- 
nian takes the form (exercise) 



H-liN = Ak ( a U a W ~ «-ki«Lki) + + 



A 2 



A 2 

= A k«L«ka + ^ (Ck - A k ) H 

kcr k ^ 

This result shows that the elementary excitations or quasi-particle states, known as "Bo- 
goluibons", created by ccj^, have a minimum energy A, the energy gap. 

To determine the ground state wavefunction one simply has to identify the state which 
is annihilated by all the quasi-particle annihilation operators a ka . This condition is met 
uniquely by the state 



|g.s.) = J^tt kT «_k|l^) ^ II ( c os6» k - sin^kC^cLk^ 



where \Q) represents the vacuum state, and 

2 sin 2 6 k = 1 - |^ 
Ak 

Physically, in the limit A — > 0, sin 2 9^ — ► — £k), and the ground state collapses to the 
filled Fermi sea with chemical potential //. As A becomes non-zero, states in the vicinity 
of the Fermi surface rearrange themselves into a bound state condensate and lower their 
energy. 



Sin 2 9 k 


Cos 







^k 

Figure 5.12: Schematic diagram showing the variation of the occupancy of the momemtum 
basis states in the ground state; sin 2 6k represents the occupancy of the k states. Note 
that the wavefunction of the ground state condensate involves the occupation of basis 
states with momentum in excess of pp = y/2m\x. 

An estimate of the corresponding ground state energy obtains 

E gA = (g.s.\H - fiN\g.s.) = (& - A k ) + — , 

k ^ 



20 



One may show that the phase of A may be chosen arbitrarily, a fact to which we will return later. 
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a result which always lowers the energy below the non-interacting g = theory (exercise). 

Finally, to determine the scale of the order parameter we have to determine self- 
consistently the order parameter, viz., 



A = 5(^(g.s.|c_ ki c kT |g.s.) = sin 6> k cos 6> k 



A 



2 ^ ( A 2 + £2)1/2 



gL 3 A r-D 
—Z—VW 



gL 3 Av(n) sinh _1 (cj D /A) 



where ^ k — > J d£ v(£) and denotes the density of states. (It is left as an exercise 
to show that a minimisation of the g.s. energy obtains the same self-consistent equation 
for the order parameter.) Here we have assumed that the pairing interaction g extends 
over an energy scale set by up- Physically, for pairing mechanisms which derive from the 
exchange of phonons, this energy scale is set by the corresponding Debye frequency, the 
maximum energy phonon. Finally rearranging this equation, one obtains 



sinh(l/ gL 3 v(fx)) 



~2ujd exp 



This completes our formal investigation of the BCS transition from the mean-field 
Hamiltonian. How does the same transition emerge from the corresponding field theory? 
In the following, we will develop a theory of superconductivity from the coherent state 
path integral for the quantum partition function. 



5.4.2 Superconductivity from the Path Integral 



To investigate the BCS transition within the framework of the coherent state path integral, 
it is convenient to abandon the long-ranged pairing Hamiltonian considered above and 
introduce a space-local attractive interaction contained within the BCS Hamiltonian, 



H 



BCS 



d d r 



5Z C -Wf^ c -W ~ 9c\(r)c\(r) Cl (r)c^r) 



Expressed in the form of the coherent state path integral, the corresponding quantum 
partition function takes the form 



D(ip, ip) exp ^ — / dr 

a.p.b.c. | Jo 



d d r 



9r + — - V ) ~ g^lpllpl^ 



where ip(r,r) represent anticommuting or Grassmann fields and the nemonic a.p.b.c. 
denotes anti-periodic boundary conditions. As usual the quartic interaction of the fields 
prevents the partition function from being evaluated explicitly. Moreover, anticipating the 
existence of a transition of the electron gas to a condensed phase in which electrons in the 
vicinity of the Fermi surface are paired, we can expect that a perturbative expansion in 
the coupling constant g will be inadequate. Motivated by the mean-field theory discussed 



Concepts in Theoretical Physics 



5.4. SUPERCONDUCTIVITY 



175 



above, we will instead introduce a bosonic field A to decouple the interaction and which 
will have the physical significance of the complex order parameter. 

The decoupling is arranged using a Hubbard-Stratonovich transformation 

'1 



o9 J ' dr j 'o CVi/> T V>iV>iV>T 



/ 



D(A*,A)exp^ - dr d 



9 



A| 2 -(A*^ t + A^) 



where A(r, r) represents a dynamically fluctuating bosonic or complex field with a sym- 
metry that reflects that of the bilinear V'j.V'T- Taking A to be homogeneous in space and 
time, the quantum Hamiltonian corresponding to the action coincides with that of the 
mean-field Hamiltonian considered in the previous section. Motivated by that analysis, 
we turn to the Nambu spinor representation 



$1 



wherein the quantum partition function assumes the form 
D(4>,iP) [ D(A*,A)exp\- / (It / d" r 



9 



lAP + ^CT 1 ^ 



[Gj p) ]- 
-A* 



-A 

where [G^] -1 = d T + (e^ — //) and [G^] -1 = d T — (ep — /j) represents the non- interacting 
Green function of the particle and hole respectively, and Q is known as the Gor'kov 
Green function. 

Being Gaussian in the fermionic fields, the functional integral over the Grassmann 
fields can be performed straightforwardly, and yields the formal expression 



= J D(A*,A)exp - J dr J 



cfV-|A| 2 + lndetCT 
9 



where we have written det Q~ x = expflndet Q~ x \. By introducing a Hubbard-Stratonovich 
decoupling of the local interaction, we have succeeded in expressing the quantum partition 
function as a path integral over an auxiliary bosonic field A. Further progress is possible 
only within some approximation. Empirically, we know that the superconducting tran- 
sition is second order, i.e. the order parameter A develops a non-zero expectation value 
below a critical temperature T c growing continuously from zero. At temperatures T <^ T c , 
spatial and temporal fluctuations around the expectation value A can be treated as small. 
In this limit, the action can be treated within a mean-field approximation where the par- 
tition function is dominated by the saddle-point configuration of A. The saddle-point 
analysis, which leads to a self-consistent equation for A known as the Gap Equation, 
is left as an exercise in Problem Set 4. Instead we will focus on temperatures T ~ T c in 
which an effective action for A can be obtained. 

In the vicinity of the transition temperature T c , the order parameter A is expected to 
be small. We are therefore at liberty to look for a perturbative expansion in powers of A. 
Setting 



Qo 1 



I -Go 



A 

A* 
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where, by definition, Q = Q(A = 0), and expanding the action to second order in A, 21 
the quantum partition function takes the form (exercise) 

Z = Z Q J D{A\A)e~ s , S = ^ [- + U(u n , q) 



5 



|A^, q | 2 + 0(A 4 ), (5.16) 



with (see Fig. 5.13a) 



n(w m ,q) = -^J^G'J p) (a; n ,k)G'J h) (c i ; n + c i ; m ,k + q). 



-k-q, -co n - co 




k,co n 

(a) 




Figure 5.13: Diagrammatic representation of (a) the response function II(u; m ,q), and (b) 
the quartic vertex. 



An instability of the electron gas towards the formation of a paired condensate is 
signalled by the appearance of a non-zero expectation value of the anomalous average 
(ijj^ijjl) (i.e. A ^ 0). Moreover, intuitively, we would expect the action to be minimised 
by a spatially and temporally uniform field configuration of A. Applying this Ansatz, we 
are led to consider a gradient expansion of the action in powers of A. 

Neglecting temporal fluctuations altogether, a gradient expansion in powers of q ob- 
tains 

K > 

U(u m , q) = 11(0, 0) + iq 2 +0(iu m , q 4 ), 

where, by symmetry, we have made use of the fact that II = II(q 2 ). Substituting this 
expansion into Eq. (5.16), and returning to the real space representation, we obtain the 
static effective action 



S[A] =(3 d d r 



^|A| 2 + ^|<9A| 2 + M |A| 4 + 



(5.17) 



where t/2 — g + 11(0,0) represents an effective 'chemical potential' for A, and u > 
represents the constant coefficient associated with the quartic vertex (see Fig. 5.13b). 
(The calculation of the coefficients K and u is left as an exercise in Matsubara frequency 
summations! However, in the following, we will not need to know their form explicitly.) 



^cre we have made use of the identity In dct Q 1 = tr In Q 1 to form the expansion 



trlnCT 1 = trlnCJo" 1 



E-tr 
n 

n=l 



Go 



A 

A* 
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Figure 5.14: Dependence of the effective action Sq on the complex mean-field order pa- 
rameter A. At t > 0, fluctuations of A are energetically unfavourable whilst for t < 0, the 
action becomes unstable towards the formation of a non-zero expectation value of A. Also 
shown is the expectation value A of the order parameter as a function of temperature. 



5.4.3 Gap Equation 

Applying a mean-field approximation, i.e. focusing on the spatially homogeneous 
component A , the effective action takes the form 



^-'lA l 2 + dA I 4 
__-|Ao| +u\A \ . 



(5.18) 



In particular, when the effective chemical potential t becomes negative, the expectation 
value of (Ao) becomes finite (see Fig. 5.14). Summing over momenta, and making use of 
the identity jd J2k = f™. d£v(£) where, as usual, v(£) denotes the density of states at 
energy £, one obtains 



n(o,o) 



-y 



-y 







E 



(3uj n 



Recalling that the attractive interaction of the electrons was mediated by the exchange 
of phonons, we note that the Matsubara summation should be cut-off at the scale of the 
Debye frequency. Setting ujd = (2n max + l)ir//3, we obtain 



n(o,o)~-vGu) E 



2n + 1 



dn 



2n + l 



Equating this result with 1/g, we deduce that the electron gas becomes unstable towards 
the formation of a pair condensate when 



T <T C = oj£> exp 



1 
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Substituting this result into the expression for t we find that, in the vicinity of T c , 



i.e. one may think of the parameter t as a 'reduced temperature'. 

Finally, taking the partition function to be dominated by the minimum of the mean- 
field action (5.18) Z ~ exp[— S (\ A|)] (i.e. applying a saddle-point approximation), we 
find a spontaneous breaking of the U (1) symmetry of the complex order parameter — i.e. 
A has a magnitude (see Fig. 5.14b), 



with arbitrary but constant phase. This situation is reminiscent of the Heisenberg fer- 
romagnet: the complex order parameter is isomorphic to a 'two- component' spin where 
the phase degree of freedom is mirrored in the orientation of the moment. According 
to the Mermin- Wagner theorem, the breaking of the continuous U(l) symmetry should 
be accompanied by the appearance of massless Goldstone modes. The stability of the 
mean-field solution towards fluctuations is governed by the effective action S[A] (5.17). 

To summarise, the quantum partition function of an electron gas subject to a local 
attractive pairing interaction has been cast in the form of a quantum field theory involving 
a complex scalar field A(r, r) whose expectation value is connected to the anomalous 
average (ip^ipi). A gradient expansion of the effective action in powers A reveals an 
instability of the electron gas towards the formation of a spatially and temporally uniform 
pair condensate. Fluctuations of the order parameter A around its mean-field expectation 
value are described by a low energy effective action (5.17). 

In fact, interpreted as an effective Free energy, this result might have been guessed 
on purely phenomenological grounds: indeed, identifying the anomalous average as an 
appropriate order parameter, A Eq. (5.17) is consistent with a gradient expansion of 
the Free energy (in powers of A) compatible with symmetry and the observed temperature 
dependences. The phenomenology of superconductivity (expressed by the free energy of 
(5.17) and known as the Ginzburg-Landau Theory) anticipated the microscopic theory 



5.4.4 f Superconductivity: Anderson-Higgs Mechanism 

> INFO. So far, our analysis of the quantum partition function associated with the BCS Hamil- 
tonian is incomplete. Indeed, no reference has yet been made to the characteristic or defining 
properties: superconductivity, and perfect diamagnetism. To discover such phenomena within 
our theory we have to generalise our approach to accommodate an electromagnetic field. At the 
microscopic level, we can do so by incorporating into the BCS Hamiltonian the canonical mo- 
mentum p — > p — eA/c, and introducing the classical action associated with the electromagnetic 
field: C = —F jXV F^ v /4, where F^ u = d fJi A u — d u A tM represents the electromagnetic field tensor. 
Leaving the formal derivation as an exercise, we can, in the spirit of Ginzburg-Landau theory, 

Concepts in Theoretical Physics 





t > 0, 
t < 0. 



of BCS 



5.4. SUPERCONDUCTIVITY 



179 



guess the appropriate form of the resulting (time-independent) action: 



S = d d r 



1|A| 2 + y \{d + i2eA/c)A\ 2 + u\A\ A + 1(0 x A) 2 



(5.19) 



This result can be inferred from gauge invariance of the order parameter under the local U{\) 



transformation tp — > e 



ie<f>(r)/c 



A(r) 



D 2ie<f>(r)/c 



A(r) — see below — (or indeed obtained from 



first principles by generalising the derivation presented above). The apparent 'doubling of the 
electric charge' can be interpreted as reflecting the effective charge of a Cooper pair. 

Gauge Invariance: In the presence of the electromagnetic field, the quantum partition 
function Z = J DA J L>[A,A*]e~ 5 exhibits a gauge invariance under the transformation (exer- 
cise) 



A' 



d</>, 



AwA' = e 2ie<l>/c A. 



That is, under such a transformation, the action remains invariant. Therefore, by gauge fixing 
d<p = A, the phase of the order parameter can be eliminated from the action and the effective 
action takes the form 

S = f3 J ' d d r 1|A| 2 + ^ (<9|A|) 2 - ^A 2 + u|A| 4 + 1(0 x A) 2 

where m 2 = 4e 2 K\A\ 2 /c 2 . As a result, we find that the massless phase degrees of freedom 
(f) have disappeared from the action! They have been subsumed into the longitudinal mode 
of the vector field A, which has itself become massive. This is an example of the celebrated 
Anderson-Higgs mechanism: below the transition temperature, the Goldstone bosons (in 
this case 4>) and the gauge field (in this case the electromagnetic field — the photon) conspire 
to create massive excitations, and the massless excitations are unobservable. 

It is instructive to interpret this result from the saddle-point or mean-field equations of 
the motion of the order parameter and vector potential. Minimising the effective action (5.19) 
with respect to spatial variations of A and A, one obtains the Gross-Pitaevskii equations 
(exercise) 



—K (d + 2ieA/c) 2 + t + 4u| A| 2 



A = 



d x (d x A) = j -m 2 A, j = 2i-K(AdA - Ad A). 



Substituting the order parameter by its homogeneous mean field value and differentiating 
the second equation, one finds that d x (d x B) = — m 2 B, from which it follows that 

(d 2 - ml) B = 0, 

where B = d x A. This result, known as the first London Equation, admits B = as the 
only constant spatially uniform solution. In a uniform superconductor the magnetic field is zero 
- the Meissner effect. At the edge of the superconductor, this equation can be integrated to 
give B ~ e~ mvX showing the field to penetrate a distance m~ l — the Penetration depth - 
into the sample. 

That these equations imply superconductivity can be inferred from the time derivative of 
the current dt} = m 2 E, where E denotes an external electric field (the second London equation). 
If a uniform field is applied for a time to, a current m 2 Eto builds up. This current remains even 
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if the electric field is subsequently switched off. This contrasts with conventional conductors 
where there is a relaxation of the current. 

Superconductivity is destroyed by a sufficiently strong magnetic field. An estimate of the 
critical field can be made by studying the Gross-Pitaevskii equation 

—K (d + 2ieA/c) 2 A(r) = tA(r) 

This equation is formally equivalent to a Schrodinger equation describing a particle of charge 
2e and mass m = 1 /2K in a uniform magnetic field. The lowest Landau level is defined by the 
condition t/K = 2eB c /c. This defines the highest field where superconductivity can occur. 

Finally, we remark that, to expel a magnetic field from a sample we require an energy 
of B 2 /2 per unit volume to resist the magnetic pressure. This must be compensated by the 
condensation energy S/j3L d . If the threshold field is smaller than the critical field B c (Type 
II) magnetic field penetrates the sample in the form of flux tubes. At low temperatures the 
latter arrange themselves in a hexagonal configuration known as an Abrikosov vortex lattice. 
Superconductors where the situation is opposite are known as Type I. 
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Figure 5.15: Bitter patter of an Abrikosov vortex lattice. 



This concludes our formal discussion of the instability of the electron gas. For the sake 
of clarity we avoided a detailed discussion of the physical manifestations of the Cooper 
instability of the electron gas — it was used here merely as a vehicle for illustrating the 
general approach of the coherent state path integral. However, before leaving this section, 
a few remarks concerning the question of universality are in order. 



5.4.5 Statistical Field Theory: Ferromagnetism Revisited 

The motivation that stands behind the phenomenology of Ginzburg-Landau theory has 
profound implications that go beyond its application to superconductivity pervading all 
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areas of physics. To illustrate the generality of the concept, let us temporarily leave 
behind superconductivity and consider the classical equilibrium statistical mechanics of a 
'one-component' or Ising ferromagnet (i.e. spin degrees of freedom can take only two 
values: Si = ±1). Our previous considerations in chapter 2 have emphasised that, when 
viewed microscopically, the development of magnetic moments on the atomic lattice sites 
of a crystal and the subsequent ordering of the moments is a complex process involving the 
cooperative behaviour of many interacting electrons. However, at first sight, this picture 
seems to be at odds with the empirical observation that thermodynamic properties of 
different macroscopic ferromagnetic systems seem to be the same — e.g. temperature 
dependence of the specific heat, susceptibility, etc. Moreover, the thermodynamic critical 
properties of completely different physical systems, such as an Ising ferromagnet and a 
liquid at its boiling point, show the same dependence on, say, temperature. What is the 
physical origin of this Universality? 

Suppose we take a ferromagnetic material and measure some of its material properties 
such as its magnetisation. Dividing the sample into two roughly equal halves, keeping the 
internal variables like temperature and magnetic field the same, the macroscopic proper- 
ties of each piece will then be the same as the whole. The same holds true if the process 
is repeated. But eventually, after many iterations, something different must happen be- 
cause we know that the magnet is made up of electrons and ions. The characteristic 
length scale at which the overall properties of the pieces begins to differ markedly from 
those of the original defines a correlation length. It is the typical length scale over 
which the fluctuations of the microscopic degrees of freedom are correlated. 

Now experience tells us that a ferromagnet may abruptly change its macroscopic be- 
haviour when the external conditions such as the temperature or magnetic field are varied. 
The points at which this happens are called critical points, and they mark a phase 
transition from one state to another. In the ferromagnet, there are essentially two ways 
in which the transition can occur (see Fig. 5.16). In the first case, the two states on 
either side of the critical point (spin up) and (spin down) coexist at the critical point. 
Such transitions, involve discontinuous behaviour of thermodynamic properties and are 
termed first-order, (c.f. melting of a three-dimensional solid) The correlation length at 
such a first-order transition is generally finite. 

In the second case, the transition is continuous, and the correlation length becomes 
effectively infinite. Fluctuations become correlated over all distances, which forces the 
whole system to be in a unique, critical, phase. The two phases on either side of the 
transition (paramagnetic and ferromagnetic) must become identical as the critical point 
is approached. Therefore, as the correlation length diverges, the magnetisation goes 
smoothly to zero. The transition is said to be second-order. 

The divergence of the correlation length in the vicinity of a second order phase transi- 
tion suggests that properties near the critical point can be accurately described within an 
effective theory involving only long-range collective fluctuations of the system. This invites 
the construction of a phenomenological Hamiltonian or Free energy which is constrained 
only by the fundamental symmetries of the system — Ginzburg-Landau theory. Al- 
though the detailed manner in which the material properties and microscopic couplings 
of the ferromagnet influence the parameters of the effective theory might be unknown, 
qualitative properties such as the scaling behaviour are completely defined. 
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Figure 5.16: Phase diagram of the Ising ferromagnet showing the average magnetisation S 
as a function of magnetic field H and Temperature T. Following trajectory 1 by changing 
the magnetic field at constant temperature T < T c , the sample undergoes a first order 
phase transition from an average 'spin-up' phase to an average 'spin-down'. By changing 
the temperature at fixed zero magnetic field, the system undergoes a second order phase 
transition at T = T c where the average magnetisation grows continuously from zero. This 
second order transition is accompanied by a spontaneous symmetry breaking in which the 
system chooses to be in either an up or down- spin phase. (Contrast this phase diagram 
with that of the liquid-gas transition — magnetisation S — > density p, and magnetic field 
H — > pressure.) The circle marks the region in the vicinity of the critical point where 
the correlation length is large as compared to the microscopic scales of the system, and 
Ginzburg-Landau theory applies. 



Following this philosophy, the Ginzburg-Landau theory of the Ising ferromagnet is 
defined by gradient expansion of the effective Free energy in powers of the order parameter, 
the local magnetisation S(r). Respecting the symmetry properties of the microscopic 
Hamiltonian (translational and rotational invariance in the spatial degrees of freedom, up- 
down or Z 2 invariance in the internal spin degrees of freedom, etc), the partition function 
of the Ising ferromagnet can be expressed in the form of a functional field integral over 
different spin configurations S(r) 



/ 



DS(v) 



-pH[S(r)] 



with the effective Ginzburg-Landau Hamiltonian or Free energy functional 



PH[S(v)} 



d d 7 



K 



-S 2 + — (dSy + uS i + --- + HS 



At the mean-field level (i.e. neglecting fluctuations of the magnetisation field) a min- 
imisation of the effective Free energy (in the absence of an external magnetic field) 



V 



uS 4 



known as the Landau Free energy, leads to the average magnetisation 

t > 0, 

y/t/Au t < 0. 
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Thus a non-zero magnetisation develops when t < identifying this parameter as the 
reduced temperature t = (T — T c )/T c . 

The Ginzburg-Landau Free energy can be compared with that obtained from the mi- 
croscopic Hamiltonian for the superconductor above (5.17). Apart from the complex 
nature of the order parameter A, the action coincides. In fact, we might very well have 
circumvented the analysis of the microscopic Hamiltonian, and written the Ginzburg- 
Landau Free energy on purely phenomenological grounds. The dependence of the param- 
eters K, u, etc., on the microscopic or material properties of the system would have been 
unavailable, but the nature of the critical point and the physical properties associated 
with the transition would have been accessible. Indeed, as mentioned above much of the 
phenomenology of conventional superconductivity was developed and understood in this 
framework even prior to the discovery of the BCS theory. 

From this result, we can draw important conclusions: Critical properties in the vicin- 
ity of a both classical and quantum second order phase transitions fall into a limited 
number of universality classes defined, not by detailed material parameters, but by the 
fundamental symmetries of the system. When we study the critical properties of the Ising 
transition in a one-component ferromagnet, we learn about the nature of the liquid-gas 
transition! Similarly, in the jargon of statistical field theory, the superconductor, with 
its complex order parameter A is in the same universality class as the two-component or 
XY Heisenberg ferromagnet. The analyses of critical properties associated with different 
universality classes is the subject of Statistical field theory. 
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